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Abstract. Matrix models of 2D quantum gravity are either exactly solvable for matter 

of central charge c < 1, or not understood. It would be useful to devise an approximate 

scheme which would be reasonable for the known cases and could be carried to the unsolved 

cases in order to achieve at least a qualitative understanding of the properties of the models. 

The double scaling limit is an indication that a change of the length scale induces a flow in 

the parameters of the theory, the size of the matrix and the coupling constants for matrix 

models, at constant long distances physics. We construct here these renormalization group 

equations at lowest orders in various cases to check that we reproduce qualitatively the 

properties of c < 1 models. 
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1. Introduction 

The matrix model representation of two-dimensional quantum gravity |l| has led to 
explicit solutions for minimal conformal fields couped to gravity || and to beautiful connec- 
tions with various branches of mathematical physics such as integrable flows or topological 
field theories. However this approach has not been of real help for understanding the diffi- 
culties which arise when the central charge c of the matter field is larger than one although 
it is very easy to write matrix models for those cases as well. Indeed if the topological 
expansion of field theories in dimension larger than one is notoriously difficult to han- 
dle, one can write (for instance) a matrix model for n Ising spins coupled to a randomly 
triangulated surface, with a central charge c = n/2, as an integral over 2 n matrices ||. 
Unfortunately these models are not solvable (up to now) for n > 1. 

Therefore we are in the unfortunate situation of either solving exactly (for c < 1) or not 
understanding at all what is happening. We are not even guaranteed that a matrix model 
candidate to describe a c > 1 model has anything to do with the continuum description, 
although we should keep in mind that (i) it is not yet clear that there is a continuum theory 
at present for c > 1, (ii) we have invested enough efforts in matrix models to try to find 
out whether they do have a continuum limit for c > 1; is there still a double scaling limit 
with some new type of exponent, or no continuum limit at all? If there is a continuum 
limit is it still a string theory in non-critical dimension? 

We have no answer to provide to these important questions, but we would like to 
devise a method which may be only approximate for the soluble cases c < 1, but which 
would allow us to understand at least qualitatively c > 1. The attempt that we present 
here is based on simple minded analogies with critical phenomena. We note first that 
the central result of matrix models is the existence of the "double scaling limit" ||, i.e. a 
continuum limit with critical exponents which describe how the two coupling constants of 
the theory have to be tuned to reach this limit. These simple scaling laws for c < 1, with 
logarithmic deviations at c = 1 [|J are reminiscent of the theory of phase transitions for 
dimensions larger or equal to four. There we know that a non-trivial IR fixed point governs 
dimensions smaller than four and that we have to rely on approximate methods such as 
thee or 1/iV expansions. Is there a similar phenomenon for matrix models? This question 
makes sense only if we understand renormalization group flows. Indeed if we had such flow 
equations, we know that the scaling laws and the exponents, which characterize the double 
scaling limit, would arise automatically. Therefore we have to understand how the two 
coupling constants of the theory, the string coupling constant and the cosmological constant 
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(mapped respectively into the size N of the matrices and the matrix coupling constant g), 
evolve under a rescaling of the regularization length introduced in the triangulation of the 
worldsheet. We thus expect that a change iV i-» N + 5N can be compensated by a change 
g i— > g + Sg with the same continuum physics. We shall see that for matrix models this 
indeed the case, at the expense of enlarging the space of coupling constants very much as 
in the Wilson's scheme of integration over the momenta in the shell A — dA < \p\ < A. 
Let us review first briefly the c < 1 scaling laws [|J: the singular part of the string 
partition function satisfies 

Z = A 2- 10 f( AN 2/-yA (-Q 

with 

A = g c - g (2) 

and 
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7i = z - 7o = — 



25-c+ v/(l-c)(25-c) . (3) 

The string susceptibility exponent at genus h is 

lh = 7o + h 71 . (4) 

The relation 

7o + 7i = 2 (5) 

independent of the explicit values of 70 and 71, is easily obtained from the consideration 
of the torus. 

Assume now that we have indeed a Callan-Symanzik like differential equation for this 
string partition function which reads 

[Nd/dN - (3(g)d/dg + >y(g)]Z(N, g) = r(g) (6) 

with a fixed point g* given by f3 (g*) = and f3' (g*) > 0. It is then elementary to verify 
that we recover the scaling law ([1]) with 

ll = HW)' 70 = 2_ ^)' <7) 

With several coupling constants gk, the scaling exponents are given by the eigenvalues of 
the matrix 

n w = H (g m m ) ■ (8) 
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Note that the scaling law (||) requires that 7 (<?*) = 2. 

These flow equations will be constructed by integrating out one line and one row of 
an (TV + 1) x (N + 1) matrix, thereby reducing it to an N x N matrix. In this process we 
shall prove that the matrix partition function Cn(9) fulfills the equation 

( N+1 (g) = [X(g)] N2 Xn {9') (9) 

with 

9' = 9 + ^M + 0(jp) (10) 

and 

It follows immediately that the string partition function 

Z(N,g) = ^ln( N (g) (12) 

satisfies (|>D with 

1(9) = 2, (13) 

and therefore the scaling law ([5]) does hold. 

The set-up of this article is the following: we first discuss at lowest order pure grav- 
ity, then multicritical points of one-matrix models, then a more appropriate saddle-point 
method and the new features which manifest themselves at higher orders. 

2. A Simple Perturbative Calculation 

We begin with a one-matrix <^> 4 model, which near its critical point describes pure 
gravity (c = 0). It consists of an integral over an iV x N hermitian matrix (p^: 

Cn(9) = / d(j) N exp [S N (0jv, g)} , (14) 



with an action 



Sn [<t>N, 9} = Ntr 



1 A? > 9 ^ 

2^JV + I^n 



(15) 



From the exact solution we know that the double scaling limit is reached in the vicinity 
of 

9c = -1/12 (16) 
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with an exponent 



7i = 5/2. 



(17) 



The matrix (pN+i is parametrized in terms of an iV x N submatrix (p^, a complex N- 
component vector v a , and a number a: 



>N+l 



<PN V a 

* 

i> a 



(18) 



but one verifies easily that all the terms involving a are of relative order 1/N and can be 
dropped in the continuum limit; in other words we can set a = 0. We then define 

exp[-S' N (c/ ) ' N ,g')} = X N (g)Jd N vd N v* exp [S N+1 (<t> N ,v,g)} , (19) 

in which <p' N is obtained after a rescaling which normalizes the coefficient of (j) N to N/2 as 
in®. 

An easy calculation yields 



S N+1 ((P N+1 ) = (N + 1) 



tT I l<t>N + l<t>i \ < ' 



+ (N+l)g 



v*<f> 2 N v + -(v*,vf 



(20) 

after which we expand the exponential to first order in v, perform the easy integrations 
over the w a 's and re-exponentiate to get a new effective action 



We then rescale 4>n to 



1 r ( ^ + 9 -cj>% 



>N = P9n 



+ gti( 



>N ■ 



with 



P=1 -^2N- + {n1) 



(21) 



(22) 



(23) 



(we are dropping from now on all terms of relative order 1/N 2 as well as higher orders in 
g) . This rescaling gives as coefficient of <p 



g' = g-^{g + ±g 2 ) 



(24) 



and for the prefactor in the partition function 



A(s) = l-i(l + 3<7). 



(25) 



Therefore we have established at first order the renormalization group equation (||) with 

(3(g) = -g-4g 2 + (g 3 ) (26a) 

7(3) = 2 (266) 

r(g) = -1 - 3a + O [g 2 ) . (26c) 

There are two fixed g* = or —1/4, but we have to select the repulsive fixed point 

9* = -1/4 (27) 

since the pure gravity exponents are obtained only when we "tune" the cosmological con- 
stant o near its critical value o c ; the true g c is —1/12 in the exact theory instead of our 
first approximation —1/4. Since at this order j3' (g c ) is equal to one, this calculation gives 

7i = 2, (28) 

to be compared to the exact value 5/2. 

We shall return to this c = case in a more elaborate calculation below, but we first 
consider within the same approximation the multicritical points of one-matrix models. 

3. Multicritical Points 

In order to study the existence of Kazakov's multicritical points J7[ within one-matrix 
models, we allow for an arbitrary polynomial in the action 

oo 

S [N, g k ]=N J2 §tr<^ (29) 

with gi — 1. 

The parametrization fll8f) of 4>n+i gives a number of terms, since 

tr <ffi +1 = tr 4$ + 2kv*(ff- 2 v + O (v 4 ) . (30) 

The crudest calculation consists in keeping these quadratic terms in v, expanding to first 
order in the Ofc's, integrating over the v's and reexponentiating. This gives, up to a nor- 
malization constant, 



Ov+i = / d(/> N exp < -(N + l)tr 



1 i g2 Ii2 , V^ f 9k 9k+i \ ,2k 



(31) 



Rescaling <f> in order to set the coefficient of 4> 2 N to AT/2, we obtain 

9k = 9k + jj [~ (2kg 2 + k-l)g k + 2kg k +i] , (32) 

and thus 

P k = - (2kg 2 + k-l)g k + 2kg k+1 . (33) 

It is easy to verify that there is a multicritical point of order m given by 



for m < k. 
The corresponding multicritical potential is 

1 - I I 




S* N = Ntr 



,2 \ m- 
2m 



(35) 



which is of course approximate; for m = 2 it gives x 2 /2 — x 4 /16, instead of the exact 
@ x 2 /2 - x 4 /48; for m = 3, x 2 /2 - x 4 /12 + x 6 /216 instead of x 2 /2 - x 4 /12 + x 6 /180. 
Although it is only approximate, we do find the right sequence of multicritical points 
with their characteristic features of being alternatively unbounded below for even m's and 
bounded for odd m's. 

The critical exponents are related to the matrix (|8|) of derivatives of the (3- functions at 
the fixed points. For the m-th multicritical points this matrix is real, (m — 1) x (m— 1) and 
non-symmetric. Remarkably enough this matrix turns out to have (m — 1) real eigenvalues 
which are 2/m, 3/m,...,m/m. For an arbitrary direction of approach of the fixed point the 
leading critical exponent is the largest one, namely one, which corresponds in fact to the 
c = exponent. If we choose a direction in the space spanned by the (m — 2) subleading 
eigenvalues we shall see in general as leading exponent (m — l)/m. In order to tune the 
m-th multicritical point we have to suppress the (m — 2) largest eigenvalues by choosing 
a direction in the g k — g\ space along the eigenvector of the smallest eigenvalue, namely 
2/m. This gives for the m-th multicritical point 

7i = m (36) 

whereas the exact value is m + 1/2. This trivial calculation is thus in good qualitative and 
semi-quantitative agreement with the exact answer. 
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4- A Better Calculation 

An attempt to expand to higher orders in the coupling constant in order to improve 
the above calculations reveals that we have not done it in a systematic way. For instance 
if we follow the same procedure to next order and try a "two- loop" calculation, we find 
that the /3-function at order g 2 receives contributions from higher orders. In fact we have 
to integrate first over the iV complex v's without perturbation theory, and then expand in 
order to reach the large N, double-scaling limit. To this effect we write the effective action 

1 



Sn+i [<(>n, v a , cr] = (N + 1 

+ (N + l)g 



tv I 2^ + f^ I -'■ '''"•'' 



* / 2 i * 

V <p N V + (TV ■ V 



a- 



(37) 



which, upon integration over a, is equivalent to (|20|). We next integrate over the u's and 
obtain 

■2 



ie& 



1 



9^ 



'• '• I 2<Pn + ^n 



.9° 



(N + l)^+trln(l + gcr + g<p 2 N ). (38) 



S%[<1> n ,<t] = (N + 1) 
In the large A^ limit, the single cr-mode is fixed to its saddle-point value given by 



o-= -^tr (1 + ga + g(p 2 N ) . 



(39) 



This gives cr as a functional of the traces of the even powers of (p, involving also non- 
linear terms as (tr<^) and so on. Therefore the original action has to be enlarged to 
accomodate such terms. However is we expand again in powers of g these new interactions 
do not yet appear at order g 2 since: 



a = 1 - (1 + 1 2 ) g + (2 + 3t 2 + U) + O (g 3 ) 



(40) 



in which 



tk = N tTi 



and thus, dropping all terms of order g s 



S 



eff 



N 



±(N + l)+g-g 2 



t 2 + 



9f N+1 \_9_ 

V } 2 



U. 



(41) 



(42) 



After a rescaling of 0jv to enforce to N/2 the coefficient of t% we find a new (p 4 coupling 
constant 

(43) 



g' = g-jj{g + Sg 2 )+o{ g *) 



i.e. 



P(g) = -g-6g 2 + 0(g 3 ). (44) 

The fixed point g* = —1/6 is a slight improvement over the previous calculation, but at 
this order the exponent 71 does not change. 

A calculation at next order requires two additional coupling constants, namely that 
■ (p® and of (tr </> 2 ) ; the renormalization flow takes pi 
The result of this calculation will be reported elsewhere. 



of tr </> 6 and of (tr (/> 2 ) ; the renormalization flow takes place in a three dimensional space. 



5. Conclusion 

Simple renormalization group transformations reproduce qualitatively and semi- 
quantitatively the results of matrix models. In principle we could apply readily the same 
technique to matrix models with c larger than one. However manifestly the method has 
first to pass a few non-trivial tests: can it be systematically improved at higher orders for 
simple one-matrix models? For c = 1 we should find that the fixed point is a double zero 
of the (3- function in order to account for the logarithmic deviations to scaling ||. The 
situation is still unclear, but we believe that it is worth exploring the possibilities of this 
approach further. 
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